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^ ■ Abstract 
^— > ■ 

With an aim to investigate the effect of externally imposed body acceleration and mag- 
'— — netic field on pulsatile flow of blood through an arterial segment having stenosis is under 
I ■ consideration in this paper. The flow of blood is presented by a unsteady micropolar 

> ; 

\ fluid and the heat transfer characteristics have been taken into account. The non-linear 

. equations that governing the flow are solved numerically using finite difference technique 

by employing a suitable coordinate transformation. The numerical results have been ob- 

OV 

served for axial and microrotation component of velocity, fluid acceleration, wall shear 
stress(WSS), flow resistance, temperature and the volumetric flow rate. It thus turns out 
■ that the rate of heat transfer increases with the increase of Hartmann number H, while 

•i-H ■ 

^ ■ the wall shear stress has a reducing effect on the Hartmann number H and an enhancing 

effect on microrotation parameter K as well as the constriction height 5. 
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1 Introduction: 

It is well known that fluid dynamical behaviour of blood through an arterial segment having 
stenosis play a vital role in cardiovascular disease. The narrowing in the artery, commonly re- 
ferred to as stenosis, is a dangerous disease and is caused due to the abnormal growth in the 
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lumen of the arterial wall. Stenosis may be formed at one or more locations of the cardiovas- 
cular system. As a result of such undesirable formation at the endothelium of the vessel wall, 
reduction of regular blood flow is likely to take place near the stenosis. If this disease takes 
a severe form, it may lead to stroke, heart attack and various cardiovascular disease. It has 
been observed that whole blood, being a suspension of blood cells in plasma, behave as a non- 
Newtonian fluid. Several attempts [1-8] have been made to understand the flow characteristic 
of blood through arteries under various assumptions. But most of these studies have failed to 
estimate on an important role in blood rheology is the motion of erythrocytes, white blood cells 
and platelets in plasma. 

Eringen [9, 10] introduced the theory of micropolar fluid, which exhibits microscopic effects aris- 
ing from the local structure and microrotation of fluid microelements are considered to be rigid. 
A subclass of microfluids such as liquid crystals, suspensions and animal blood are treated as mi- 
cropolar fluid, which takes into account the rotation of fluid particles by means of an independent 
kinematic vector called the microrotation vector. With this end in view, several investigators 
[12-14] considered micropolar fluid as a blood flowing through an arterial stenosis. They sug- 
gested that blood exhibits remarkable non-Newtonian properties, since blood is a suspension 
of neutrally bouyant deformable particles in viscous fluid. Devanathan and Parvathamma [15] 
developed a mathematical model for steady flow of micropolar fluid through an arterial stenosis 
and their results indicated that the axial velocity is highly affected owing to the presence of 
microelements in the fluid. However, all the above studies did not considered the effect of body 
acceleration. Although, under normal physiological conditions, flow properties of blood depends 
upon pumping action of the heart, which in turn produces a pulsatile pressure gradient through 
the cardiovascular system. Although, there may be some another important situation under 
which human beings experiences whole body acceleration. It may occur for a short or a long 
period of time. Such situations arise day to day life of human being e.g., driving a vehicle, 
flying in an aircraft, due to which there may occur, the serious health problems like headache, 
loss of vision, increase of pulse rate and hemarage in face, neck and eye socket. Keeping this 
in mind, several investigators [16-24] carried out their studies pertaining to blood flow through 
stenosed arteries with a periodic body acceleration. Chakravarty and Sanniggrahi [25] studied 
theoretically the flow characteristics of blood through an artery in the presence of multi-stenosis 
when it is subjected to whole body acceleration. Wherein they considered a flexible cylindrical 
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tube containing a homogeneous Newtonian fluid with variable viscosity characterising blood. 
But all the above studies neither consider the effect of body acceleration on pulsating blood flow 
nor includes the effect of magnetic field when blood is considered to be micropolar fluid model. 
Since blood consists of suspension of red cells containing hemoglobin, which contains iron oxide, 
it is quite apparent that blood is electrically conducting and exhibits magnetohydrodynamic flow 
characteristics. If a magnetic field is applied to a moving and electrically conducting fluid, it will 
induce electric as well as magnetic fields. The interaction of these fields produce a body force 
known as Lorentz force, which has a tendency to slow down the fluid motion. Such analysis may 
be useful for pumping of blood and magnetic resonance imaging (MRI). The magnetic field can 
also be used for controlling blood flow during surgery. Many authors [26-29] have investigated 
the flow of blood through arteries in the presence of magnetic field under different situations. 
The steady MHD flow of a viscous fluid in a slowly varying channel in the presence of a uniform 
magnetic field was explored by Misra et al [30]. Recently Misra and Shit [31] investigated the 
flow and heat transfer of a viscoelastic electrically conducting fluid under the action of a mag- 
netic field. They made an observation that the temperature of blood increases as the strength 
of magnetic field increases is likely to be useful in the development of new heating methods. 
With the above discussion in mind, a good attempt is made in the present theoretical investiga- 
tion to examine the pulsatile blood flow characteristics through a stenosed artery experiencing 
periodic body acceleration by treating micropolar fluid model for blood in the presence of mag- 
netic field. The motivation of this study is to analyze the effect of microrotation component and 
axial velocity when an external magnetic field and body accelerations are applied. In addition, 
in view of the information stated above, we have also performed a heat transfer phenomena of 
the problem in question. 

2 Mathematical Formulation 

Let us consider a fully developed unsteady, laminar, incompressible and axially symmetric two 
dimensional flow through a stenosed artery. An external magnetic field is applied to a pulsatile 
flow of blood, which is characterized by a micropolar fluid model under periodic body accelera- 
tion. Let (r*,6*,z*) be the co-ordinates of the material point in a cylindrical polar co-ordinate 
system, of which z* is taken as central axis of the artery. We denote u* and v* the velocity 
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components of blood along the axial and radial directions respectively and w* the component of 
microrotation. For axisymmetry, we assumed that the flow variation is independent of 9*. Ow- 
ing to the above considerations, the governing equations of the present problem may be written 
in the form, 

du* v* dv* 

du* ^du* *du* dp* k)(^ 2u * 1 d 2u *\ U® w * w * 

-aB 2 u* + pG*(t*), (2) 

dv* ^dv* *dv* dp* d 2 i;* 1 V * d 2 v* ^dw* 

* Qz * ^ * ^ ^ f * df ^ t* * ^ Qz * ^ Q z * ' 



,<9w* ^dw* *<9-uA , , * (9-u* <9-y* N .(9 2 tu* 1 <9u>* u>* <9 2 w*, 
l_ M * |_ ^* j — k(2w -\ ) + t( I I 

1 ^ Q z ^ Qt* ^ df * d z * * ?^ ^ Qt' ^ ^ * 2 Q z ^ ^ ' 



<9t* cfc* <9r* pCp <9r* 2 r* <9r* <9z* 2 pC p 

where p represents the blood density, p the blood pressure, \i the viscosity of blood, k the 
rotational viscosity, a the electrical conductivity, B the applied magnetic field intensity, 7 the 
spin gradient viscosity is assumed to be 7 = j(fi + k/2), j the micro-gyration constant, T* is the 
temperature, C p the specific heat and kq denotes thermal conductivity. In the above, we have 
neglected the effect of induced magnetic field due to low magnetic Reynolds number (R m << 1). 
For t* > 0, the flow is assumed to have a periodic body acceleration G*(t*) appeared in equation 
(2) has the expression of the form 

G*(t*) =a*(cosuj b t* + <f) g ) (6) 

where a* denote the amplitude of body acceleration, u>j, the frequency of body acceleration and 
(j) g the phase difference. 

Let us introduce the following non-dimensional variables to put the above equations in dimen- 
sionless form: 



Using these dimensionless variables the equations (l)-(5) becomes 

® u _|_ v _|_ ^ v o (7) 
dz r dr 

2 ,du du du, dp T ^,d 2 u ldu d 2 u, T ,,dw w, tt2 , , 

2 dv dv dv, dp _ K)(^ 2v V ® 2y \ ^ w (g\ 

dt dz dr dr dr 2 r dr r 2 dz 2 dz 

2 T( dw dw dw du dv d 2 w ldw w d 2 w 

dt dz dr dr dz dr 2 r dr r 2 dz 2 



m_ m_ de__ i a 2 e im d^e_ e c h 2 2 (1V 

dt dz dr a 2 P r dr 2 r dr dz 2 a 2 ' 



where the non-dimensional parameters appeared in equations (8)-(ll) are defined as the ratio 
of viscosity K — ^, the material constant m = Womersley number a = Rq^J^-, Hartmann 
number H = RnB .f^, Prandtl number P r = and Eckert number E c = n !t, R % , . 
Using the non-dimensional quantities ao = pR ^ , b — ^, expression for the body acceleration 
takes the form 

G(t) = a cos(bt + <p g ) with t > (12) 

Also the non-dimensional form of pressure gradient and wall motion are assumed as 

dP{t) 



dz 

and 



K + K p cos{t) (13) 



R(t) = R[l + K r sin(t + r )] (14) 

where R, K r , r , K, and K p stand for the mean radius, amplitude of arterial wall motion, 
the phase difference, constant amplitude of pressure gradient and amplitude of the pulsatile 
component gives rise to systolic and diastolic pressure. It is clear from equations (8)-(ll) that 
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when K and m are set to zero, the system of equations reduces to a classical Newtonian fluid 
model. 

Let us simplify the equations of motion by using long wavelength approximation << 1), 
where the lumen of the arterial radius R is sufficiently smaller than the wavelength A of pressure 
wave. Under this assumption the axial viscous transport terms are negligible from (8)-(ll). 
Therefore the equation (8) reduces to simply |^ = and others becomes 

® u _|_ v _|_ ® v — o (15) 

dz r dr 

9 ,du du du. dp ,„ Tr .,d 2 u ldu, T ^,dw w. tt0 „. . 

" < m + "97 + "*) = -£ + ^ + A + ;<F> + ~dr + 7> - + G « < 16) 
° J( « + u & + " * » = " A (2w + * - al> + ro( a? + ; * " (17) 

<E + U ^L + V ?1- 1 (& 6 | | EcR2 (18) 

<9t cte <9r a 2 P r dr 2 r dr a 2 

The boundary conditions for the present problem are described as follows: 

Along the central line of the tube (at r = 0) the radial velocity, axial velocity gradient and spin 

of microrotation vanishes, 

du 89 

i.e, at r — 0, v — w — — = — = 0. (19) 

ar ar 

and at the wall 

dR 

i.e, at r — R(z,t), u = w = = 1.0 and v = (20) 

The geometry of the stenosis is described mathematically in non-dimensional form as (cf. Fig. 
1) 



R(z,t) = R 



(5 r , /27T ((),,, 

2R { +C ° s{ ~k {z ~ 2 ))} 



(1 + K r sin(t + (f> r )); when d<z<d + l , 



= R [1 + K r sin(t + 4>r)] ; elswhere (21) 

where .R denotes the mean radius of the artery, 5 the depth of the stenosis, R is the radius of 
the artery at normal pathological state. 
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3 Method of Solution 



Since the governing equations are couple and non-linear in nature, finding analytic solution is im- 
possible. To simplify the governing equations, we introduced a radial coordinate transformation 
given by 

^ R(z,t) 

Using this transformation, the equations (15)-(18) can be written in the following form 

du £ du dR v 1 dv 



dz R di dz £R R d£ 



0, 



(22) 



<-, du 

du _qi 



dt 



R 



'.dR dR, 



U 



du 

dz 



a 2 R 2 



d 2 u 1 du 



K dw w 

+ ^R^~bJ + J' 



+- 



1 



G(t) - H 2 u - 



dp 

dz 



(23) 



dw 
= J%- 

dt R 



dw 



'.dR dR, 



dw K dv , 



K 



m 



a 2 JR 
d 2 w 



du 



.dv dR 

dt; dz 



a 2 JR 2 



de 



1 dw w 



(24) 



dt R 



'.dR dR, 



de 

-u— + 



1 



dz a 2 P r R 2 



d 2 e \de_ 
W 2 + Wi 



+ 



E r H< 



The boundary conditions (19) and (20) are also transformed to 

du d9 . 
v = «, = - = - = at ^ = 0. 



and 



(25) 



(26) 



dR 

u = w = 0, v — — — and 9 = 1.0 at £ = 1. 
dt 



(27) 



Let us multiply the equation (22) by £R and integrating with respect to £ between the limits 
and £, we get 

x dR R [Sdu „ 2dR /* „ 



(28) 
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Using the boundary condition (27) at £ = 1 the equation (28) yeilds 



1 dR 



2 dR r 1 



R dt 



R dz Jo 



o^T = -7^ 



1 ^du 



dz 



and has the form 



jo jo 



1 9i2 2 

R~dt n ^ } ~R~dz~ u 



(29) 



Let us choose /(£) = -4(£ 2 - 1) such that /„* = 1- 

Now equating the integrands from both sides of (29) we have 

du _ 4(£ 2 -l)dR 2 dR 



— ^u, 



dz R dt Rdz ' 

which on substitution in (28), we obtain the radial velocity in the following form 



(30) 



dR tn _ 2 .dR 



(31) 



In the next section, we present computational scheme for computing numerical solution of the 
velocity components. 



4 Computational Scheme 

Since our Physical domain is highly complex to reduce it to a simple rectangular domain, we used 
the co-ordinate transformation, presented in the previous section. To solve the equations (22)- 
(25) and (30) using a finite difference technique, we subdivided the rectangular domain into 
a network by drawing straight lines parallel to the co-ordinate axes. The solutions u(z,£,t), 
v(z,£,t), w(z,£,t) and 9(z,£,t) at any mesh point in a computational domain are denoted by 
ufj, vfj, wfj and 9fj in which 
Zi = iAz ; i = 0, 1, • • • , M; 
./Ai : J = 0,l,---,iV; 
t k = (k-l)At; A; = 1,2, 3 

where M and N are the maximum number of mesh points in the z and ^-directions respectively. 
Let V*j stand for the variables u, v, w and 9, then the spatial and time derivatives are replaced 
by the central difference as 
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9V\k _ V&+1-V&-1 (dV^k _ Vj+ij-VjLij (d_V\k _ ^,+1-2^+^-1 (9V\k _ KT- V iJ 
)i,3 2A£ ' v dz 2Az ' v <9£ 2 >i,3 (A^) 2 ' v dt >i,3 At 



v d£ >i,j 

so on 



Using these finite differences, the discretized equations for (23)-(25) and (30) becomes 



and 



u 



k+l 



u\ 3 + At[- 



- U i,j+1 U i,j-1 



'OR' 



dR k ' 



dz 



k U i+l,j U i-l,j 



i,3 / h3 



2Az 



+ 



1 + K ( u k ij+l - 2u k hj + Uj j_i | 1 u k j+1 - u k jA | K w^-wfj^ | 



(AO 5 



2A£ 



a 2 ^ v 2A£ 



(32) 
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<9z 
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2A^ 



a 2 ! ^ 



k k 

v i+ij - v i-i,j 



k„,k 



2Az 



K ( u k j+1 - u^._i + ^ dR k v iJ+1 - v_ M _ 1 | + 



a 2 R k J 



m 



2A£ 



dz i 2Ai 



'w k j+l - 2wf d + wf^ 1 w* J+1 - w\ 4 _ x 



w 



a 2 (R k ) 2 J 



+ 



1,3 



2A£ 



}, (33) 



nk ok 

est 1 "^r 1 



' (dR k k 8R ks 



fc 



\ nk _ nk 

_.,k \ _ k 

'•' 2A; 



a 2 P r {R k ) 2 



(AO 5 



2A£ 



cr 



(34) 



dz i 



dR k 



(35) 



where and (j^f) are computed at the arterial wall. 

After having determined the values of the axial velocity, radial velocity, microrotation component 
and temperature one can easily obtain the fluid acceleration F, the volumetric flow rate Q, flow 
resistance A, wall shear stress t w and the Nusselt number Nu from the following relations, 



, , i u k+1 -u k 
k+l _ u i,j u i,j _,_ K 

At " r a i,3 2Az 



_p fc + i - Z.J. ILL ^ k U i+l,3~ U i-l;3 



■I -J 



(dR k , ,.fc M^.j + l-M^,! 

R fe \dti ^ U i,j dz i J 2A£ 



Q l = 2n(R k ) 2 Jo 1 0^/0 



Nu 



k •) 



J i,N "i,N-l 



A? 
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5 Computational Results and Discussion 



The objective of the present study has been to investigate the flow characteristics of blood 
through stenosed artery under the action of a periodic body acceleration as well as an external 
magnetic field on the axial velocity by taking into account the effect of microrotation of blood 
cells. A specific numerical illustration has been taken to examining the applicability of the phys- 
iological data available in the existing scientific literatures [3,4,19-21,28,29]. The computational 
work has been carried out by using the following data: 

L = 5.0, d = 2.0, l = 1.0, 5 = 0.1, 0.25, 0.50; a = 0.0, 1.0, 2.0, 3.0; b = 1.0, <p g = 0.0, <p r = 0.0, 
R = 1.0, R= 1.0, K p = 1.46, K = 7.30, K r = 0.05, f p = 1.2, K = 0.0, 0.1, 0.2, 0.3; 
J = 0.1, m = 0.1,0.01,0.001; a = 3.0, H = 0.0, 1.0, 2.0, 3.0, 4.0; E c = 0.0002, 
and for a human body temperature, T = 310K the value of P r = 21 is considered for blood. For 
the sake of comparison, we have also examined the cases where P r = 7, 14. We used A£ = 0.025, 
Az = 0.05, At = 0.001 through out the computation. We observed that further reduction in 
A£ and Az does not bring about any significant change, which leads to the stability of the 
numerical methods. 

Figs. 2(a)-(c) give the variation of axial velocity u along the radial direction £. In Fig. 2(a) 
dotted lines indicate the axial velocity distribution for Newtonian fluid model and solid lines 
represent that of the micropolar fluid model. From this figure, we observe that the axial velocity 
decreases with the increase of the Hartmann number H . The axial velocity occurs maximum at 
the central line of the artery in all four cases. One can note from Fig 2(a) that for a sufficient 
strength of magnetic field, the profiles for both the Newtonian and non-Newtonian models are 
same. The axial velocity is found to be greater in the case of Newtonian model than that of 
the micropolar fluid model. This happens due to the presence of microrotation of blood cells, 
creates an additional viscosity called rotational viscosity, which in turn diminishes axial velocity. 
Fig. 2(b) depicts the variation of axial velocity for different values of K, the ratio of viscosity -. 
When the values of K increases, the axial velocity decreases. The results, thus obtained due to 
the increase of the microrotations of blood cells in the lumen of the artery. The axial velocity is 
minimum in the vicinity of the arterial walls, which gives rise to conclude that the microtation 
of blood cells occur in the plasma layer. The variation of axial velocity for different amplitudes 
of body acceleration is shown in Fig. 2(c). It reveals that the axial velocity increases as the am- 
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plitude of body acceleration increases. Thus in the presence of vibration environmental system 
may produce the increase in blood velocity. 

Figs. 3(a)-(c) illustrate the distribution of micro-component of velocity w along the radial 
direction of the artery. Fig. 3(a) shows that the microrotation component increases with the 
increase of the ratio of the viscosity K, which in turn produces decrease of axial velocity. It has 
been shown in Fig. 3(b) that the microcomponent of velocity w increases with the decrease in 
micropolar material parameter m. The result lies in the fact that, when m decreases, the fluid 
viscosity /x increases, thereby the microcomponent of velocity increases. One can note that the 
effect of microrotation observed significantly near the arterial wall. The variation of microrota- 
tion component is found to increase with the increase of the blockage of an artery presented in 
Fig. 3(c). In all these cases 5 = 0.1, 0.25, and 0.5 have been considered for computation and 
which are respectively represent the 10%, 25% and 50% of the blockage of the arterial segment. 

In Figs. 4(a) and (b) we present the variation of dimensionless temperature distribution at 
the throat of the stenosis and the rate of heat transfer at the arterial wall along with the axial 
direction. It is observed from Fig. 4(a) that the temperature decreases with the increase of the 
Prandtl number P r . However, the temperature increases with the increase of Hartmann number 
H. It is interesting to mention here that the temperature rapidly increases with the increase of 
H up to 4.0 (T) and beyond which, no significant change is noticed. The rate of heat transfer 
(Nu) increases with the increase of prandtl number P r , shown in Fig. 4(b). The rate of heat 
transfer is maximum at the throat of the stenosis and is minimum at the down stream of the 
stenosis than that of upstream. Thus the temperature can be increased whenever necessary by 
the application of the certain magnetic field strength. It is well known that the objective of the 
hyperthermia in cancer therapy is to raise the temperature above a therapeutic value 42° C, 
while maintaining the surrounding temperature at sublethal temperature [32]. 

The important characteristics of blood flow are the wall shear stress r w and the volumetric 
flow rate Q. Figs. 5(a)-(d) give the variation of wall shear stress r w with axial distance z and 
Time T = t/2irf p for different height of the stenosis S, different values of Hartmann number H, 
amplitude of body acceleration ao and the material parameter K. In Fig. 5(a), we observe that 
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the wall shear stress increases as the height of the stenosis increases. It is interesting to note 
from this figure that the wall shear stress significantly changes at the down stream of the steno- 
sis, whereas, at the upstream of the stenosis no change is occur. Fig. 5(b) presents the variation 
of wall shear stress r w with time T for different values of the Hartmann number H. It shows 
that the wall shear stress varies periodically with time T and the peak value of each oscillation 
decreases with the increase of H . Fig. 5(c) shows that the amplitude of the oscillation of wall 
shear stress increases with the increase of the amplitude of the body acceleration. However, 
there is no significant change on the wall shear stress for the variation of material parameter K 
(cf. Fig. 5(d)). This is due to the fact that the microelements are unable to rotate close to the 
arterial wall. Although the peak value of oscillation increases with the increase of K, but the 
change is insignificant. 

The variation of volumetric flow rate Q with time T presented through the Figs. 6(a)-(c). 
We observe that the flow rate Q decreases with the increase of Hartmann number H as well as 
the increase of the material parameter K, while the trend is reversed in the case of the amplitude 
of body acceleration G(t). In all these three cases, the flow rate Q varies periodically with time 
T. It thus turns out that under the action of a magnetic field, volume of blood flow can be 
controlled during surgeries. The volumetric flow has an enhancing effect with its peak value of 
oscillation on the vibration environmental system. 

Figs. 7(a) and (b) depict the variation of fluid acceleration F along the radial coordinate £ 
for different values of the amplitude of body acceleration and different magnetic field intensity 
given by H . For each of these cases, the fluid acceleration F is maximum at the central line and 
then it is gradually diminishes and ultimately vanishes at the wall of the artery. It is also ob- 
serve from Fig. 7(a) that the magnitude of the fluid acceleration F increases as the amplitude of 
body acceleration increases. Thus in the presence of an environmental system fluid acceleration 
enhances. However, the magnitude of the fluid acceleration F decreases with the increase of the 
Hartmann number H, shown in Fig. 7(b). One can mention here that the rate of increase or 
decrease of fluid acceleration F significantly affected by the application of magnetic field more 
than in the presence of body acceleration. 
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The variation of flow resistance A along with the height of the stenosis S has been shown in 
Figs. 8(a)-(c). It has been observed from Fig. 8(a) that the resistance to flow decreases with an 
increase in amplitude ao of the body acceleration. On the other hand flow resistance increases 
with the increase of the values of H and K, shown in Figs. 8(b) and (c). It is worth mentioning 
that in all these figures, the flow resistance A gradually increases as the depth of the stenosis 
increases. It is learn from these figures that the flow resistance remains fixed upto a certain 
height of the stenosis that is upto 25% of the blockage of an artery (called mild stenosis) beyond 
which it increases significantly. 



6 Concluding Remarks 

The present mathematical analysis is motivated towards the modelling of pulsatile blood flow 
through a stenosed artery under the combined influence of an external magnetic field and peri- 
odic body acceleration. The mathematical model is developed towards the application in blood 
flow in order to consider the effect of micro-rotation of micro-particles suspended in plasma. 
When a magnetic field is applied to the blood flow, the erythrocytes orient their disk plane 
parallel to the axis of magnetic field, creates an additional viscosity, which in turn produces 
decrease of blood velocity. This interesting observation can be seen in Fig. 2(b). The temper- 
ature of blood can be increased by the application of magnetic field. It is also concluded that 
after certain magnetic field strength, the temperature is constant. This important result can 
be useful in the therapeutic treatment of patient [32]. The present study also enables to have 
an estimate of the effects of microcomponent of velocity under the influence of periodic body 
acceleration. However, the consideration of pulsatile flow on the cases of the severe stenosis are 
the further scope of the study. 
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Fig. 1 A physical sketch of the model artery with stenosis. 
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Fig. 2 Variation of axial velocity u at the throat of the stenosis for (a) different values of the 
Hartmann number H, when K = 0.1, a = 1.0; (b) different values of K when H = 2.0, 
a = 1.0; (c) different amplitude of body acceleration a when H = 2.0 and K = 0.1. 
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Fig. 3 Variation of microrotation component w at the throat of the stenosis (a) for different 
values of material parameter K with m = 0.1, 5 = 0.25; (b) for different values of m, with 
K = 0.1, 5 = 0.25; (c) for different values of 5, with K = 0.1, m = 0.1. 
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Fig. 4 (a) Distribution of temperature 9 for different values of H and Pr. (b) Variation of 
Nusselt number Nu along the axis of the artery for different values of Prandlt number Pr. 
With a = 1.0, b = 1.0, 5 = 0.25, <p g = 0.0, r = 0.0, K = 0.1, m = 0.1, a = 3.0, H = 2.0, 
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Fig. 5 Variation of wall shear stress t w (a) along the axis of the artery for different height of 
the stenosis 5, with H = 2.0, a = 1.0, K = 0.1; (b) with time T for different values of H, 
when 5 = 0.25, a = 1.0, K = 0.1; (c) with time T for different amplitude of body acceleration 
a , when 5 = 0.25, H = 2.0, K = 0.1; (d) with time T for different values of K, when 5 = 0.25, 

a = 1.0, H = 2.0. 
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Figs. 6 Variation of volumetric flow rate Q with time T(a) for different values of H, when 
K = 0.1, a = 1.0; (b) for different values of H, when a = 1.0, H = 2.0 and (c) for different 

values of a , when K — 0.1, = 2.0. 
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Fig. 7 Variation of fluid acceleration F along the radial direction for (a) different amplitude of 
body acceleration a when H = 2.0 and (b) for different values of Hartmann number H when 
a = 1.0. 
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Fig. 8 Variation of flow resistance A with depth of the stenosis (a) for different values of do, 
when H = 2.0, K = 0.1, (b) for different values of Hartmann number H, when ao = 1.0, 
K = 0.1 and (c) for different values of K, when H = 2.0, a = 1.0. 



24 



